Chern semimetal and Quantized Anomalous Hall Effect in HgCr 2 Se 4 
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In 3D momentum space, a topological phase boundary separating the Chern insulating layers from 
normal insulating layers may exist, where the gap must be closed, resulting in a "Chern semimetal" 
state with topologically unavoidable band crossings at the fermi level. This state is a condensed- 
matter realization of Weyl fermions in (3+1) D, and should exhibit remarkable features, such as 
magnetic monopoles and fermi arcs. Here we predict, based on first-principles calculations, that 
such a novel quantum state can be realized in a known ferromagnetic compound HgCr2Se4, with a 
single pair of Weyl fermions separated in momentum space. The quantum Hall effect without an 
external magnetic field can be achieved in its quantum-well structure. 

PACS numbers: 71.20.-b, 73.20.-r, 73.43.-f 



Under broken time reversal symmetry, the topological 
phases of two-dimensional (2D) insulators can be charac- 
terized by an integer invariant, called Chern number [lj, 
which is also known as the TKNN number 0] or the num- 
ber of chiral edge states [3] in the context of the quantum 
Hall effect. 2D insulators can thus be classified as nor- 
mal insulators or Chern insulators depending on whether 
or not the Chern number vanishes. Since the Chern in- 
variant is defined only for 2D insulators, it is natural to 
ask what is its analog in 3D? Starting from a 2D Chern 
insulating plane (say at k z — 0), and considering its evo- 
lution as a function of k Zl generally two situations may 
happen. If the dispersion along k z is weak, such that 
the Chern number remains unchanged, the system can 
be viewed as the simple stacking of 2D Chern insulating 
layers along the z direction. Such 3D Chern insulators 
are trivial generalization of the Chern number to 3D, 
which is quite similar to the weak topological insulators 
in systems with time reversal symmetry. However, if the 
dispersion along k z is strong, such that Chern number 
changes as the function of k z , the system will be in a non- 
trivial semimetal state with "topologically unavoidable" 
band crossings located at the phase boundary separat- 
ing the insulating layers in k space with different Chern 
numbers 0, H| . This is due to the fact that the change of 
Chern number corresponds to a topological phase transi- 
tion, which can happen only if the gap is closed. From the 
Kohn-Luttinger theorem, we can always expect that the 
band crossings appear at the fermi level at stoichiometry. 

This Chern semimetal state, if found to exist, can be 
regarded as a condensed matter realization of (3+l)D 
chiral fermions (or called Weyl fermions) in the relativis- 
tic quantum field theory, where the field can be described 
by the 2-componcnt Weyl spinors [6] (either left- or right- 
handed), which are half of the Dirac spinors and must ap- 
pear in pairs. The band-crossing points or Weyl nodes are 
topological objects in the following senses. First, since 
no mass is allowed in 2x2 Hamiltonian, the Weyl nodes 
should be locally stable and can only be removed when a 



pair of Weyl nodes meet together in the k space. Second 
the Weyl nodes are "topological defects" of the gauge 
field associated with the Berry's curvature in momen- 
tum space 0-01 . The gauge field around the neighbor- 
hood of the Weyl node must be singularly enhanced, and 
behaves like magnetic field originating from a magnetic 
monopole [|| . The physical consequence of such a gauge 
field has been discussed in the context of anomalous Hall 
effect 0, H, E3 observed in ferromagnetic (FM) metals, 
where the Weyl nodes, if any, are always submerged by 
the complicated band structures. 

In this Chern semimetal state, we may also expect un- 
usual features like non-closed fermi surfaces (fermi arcs) 
on the side surfaces. The possible fermi arcs have been 
recently discussed from a view point of accidental degen- 
eracy, and prospected for non-collinear antiferromagnetic 
pyrochlore iridates [llj by the fine-tuning of electron cor- 
relation U. Since the correct U and the real magnetic or- 
dering is still unknown, we have to wait for its material 
realization. In this paper, we will show that such a novel 
Chern semimetal state is actually realized as the ground 
state of a known FM material HgCr2Se4, with only a sin- 
gle pair of Weyl nodes separated in the momentum space. 
We further find that the long-pursuing quantized anoma- 
lous Hall effect (QAHE) [12j4l5j, i.e., the quantum Hall 
effect without an external magnetic field, can be achieved 
in the quantum- well structure of HgCr2Sc4. 

HgCr 2 Se 4 is a FM spinel exhibiting large coupling ef- 
fects between electronic and magnetic properties 161. It 
shows novel properties like giant magnetoresistance [l 7j j , 
anomalous Hall effect [l8j], and red shift of optical ab- 
sorption edge 19]. Its Curie temperature T c is high 
(around 106~120 K), and its saturated moment is 5.64 



/is/f.u. |2Ct |21| . approaching the atomic value expected 
for high-spin Cr 3+ . Its transport behavior is different 
from other FM chalcogenide spinels, like CdCr2Se4 and 
CdCi'2S4 which are clearly semiconducting. HgCr2Se4 
exhibits semiconducting character in the paramagnetic 
state but metallic in the low temperature FM phase 17, 
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FIG. 1: (Color online) Electronic Structures of 
HgCr2Se4. (a) The total and partial density of states (DOS); 
(b) The band structures without SOC (showing the up and 
down spin parts separately); (c) The schematic understand- 
ing for the band-inversion, where the \S) state is lower than 
the |P) states in the down spin channel; (d) The band struc- 
ture after including SOC (with majority spin aligning to the 
(001) direction). The low energy states at V are indicated as 
explained in the main text. 



22| . |23[ . The spinel structure (space group Fd3m) can be 



related to the zinc-blende and diamond structures in the 
following way. If we treat the Cr 2 Se4 cluster as a single 
pseudo-atom (called X) located at its mass-center, then 
the HgX sublattice forms a zinc blende structure. There 
are two HgX sublattices in each unit cell, and they are 
connected by the inversion symmetry similar to the two 
atoms in the diamond structure. The pseudo atom X is 
actually a small cube formed by Cr and Se atoms located 
at the cube corners. The cubes are connected by corner 
sharing the Cr atoms. As the result, each Cr atom is 
octahedrally coordinated by 6 nearest Se atoms. 

Our first-principles calculations 24] confirm that FM 
solution is considerably (2.8eV/f.u.) more stable than 
non-magnetic solution, and the calculated moment (6. 
/is/f.u) is in good agreement with experiments [20I l2lj]. 
The electronic structures shown in Fig. 1(a) and (b) sug- 
gest that the system can be approximately character- 
ized as a "zero-gap half-metal" in the case without spin- 
orbit coupling (SOC). It is almost a half- metal because 
of the presence of a gap in the up-spin channel just 
above the Fermi level; it is nearly zero-gapped because of 
the band-touching around the T point just below Fermi 
level in the down spin channel. The 3d-states of Cr 3+ 
are strongly spin-polarized, resulting in the configuration 
t^e^t^ge^. The octahedral crystal field surrounding Cr 
atoms is strong and opens up a band gap between the 



hybridization with Cr-3d states, the Se-4p are slightly 
spin-polarized but with opposite moment (about -0.08 
/is/Se). The zero-gap behavior in down spin channel is 
the most important character here (Fig. 1(b)), because it 
suggests the inverted band structure around the T point, 
similar to the case of HgSe or HgTe [25|, [26[ . 

The four low energy states (8 after considering spin) 
at the T point can be identified as \P X ), \P y ), \Pz), and 
\S), which are linear combinations of atomic orbitals [27j | . 
Considering these 4 states as basis, we now recover the 
same situation as HgSe or HgTe, and the only difference 
is the presence of exchange splitting in our case. Here the 
band inversion (see Fig. 1(c), \S,l) is lower than \P, J.)) is 
due to the following two factors. First, the Hg-5d states 
are very shallow (located at about -7.0 eV, Fig. 1(a)) and 
its hybridization with Se-4p states will push the anti- 
bonding Se-Ap states higher, similar to HgSe. In addi- 
tion to that, the hybridization between unoccupied Cr- 
3d} and Hg-6s^ states in the down spin channel will push 
the Hg-6s^ state lower in energy (Fig. 1(b) and (c)). As 
the results, the \S,\.) is about 0.4eV lower than \P,i) 
states, and it is further enhanced to be 0.55eV in the 
presence of SOC. We have to be aware of the correla- 
tion effect beyond GGA, because the higher the Cv-3d^ 
states, the weaker the hybridization with Hg-Gs^. It has 
been shown that semiconducting CdCr2S4 and CdCr2Se4 
can be well described by the LDA+U calculations with 
effective U around 3.0 eV 28, 2t|. We have performed 



t^ g and manifolds. The Se-4p states (located from 
about -6 eV to eV) are almost fully occupied and con- 
tribute to the valence band top dominantly. Due to the 



the same LDA+J7 calculations for HgCr2Se4 and found 
that the band-inversion remains unless the U is unrea- 
sonably large (> 8.0 eV). The experimental observations 
of metallic behavior at low temperature for all kinds of 
samples 17, 22, 23| are strong supports to our conclusion 
for the inverted band structure. 

In the presence of SOC, the new low energy eigen states 
at T are given as ||,±|), || ,±5), If ,±5), and \S, ±i), 
which can be constructed from the |P) and IS*) states [3l| . 
similar to HgSe again. Now because of the exchange 
splitting in our case, the eight states at T are all energet- 
ically separated, with the ||,|) having the highest en- 
ergy, and the \S, — |) being the lowest. Due to the band- 
inversion, several band-crossings are observed as shown 
in the band structure (Fig. 1(d)). Among them, however, 
only two kinds of band-crossings (called A and B) are 
important for the states very close to the Fermi level. 
The crossing-A gives two points located at k z = ±k z 
along the T — Z line, and the trajectory of crossing-B is 
a closed loop surrounding T point in the k z =0 plane, as 
schematically shown in Fig. 2(a). For the 2D planes with 
fixed-fc z (k z ^ and k z ^ ±&£), the band structures are 
all gapped (in the sense that we can define a curved fermi 
level). We can therefore evaluate its Chern number C for 
each k z -fixed plane. We found that C = for the planes 
with k z < ~k z or k z > k z , while C = 2 for the planes 
with — k c z < k z <k c z and k z / 0. We therefore conclude 
that the crossing-A located at the phase boundary be- 
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FIG. 2: (Color online) Weyl nodes and gauge flux in 
HgCr2Se4 (a) The band-crossing points in the fc-space; (b) 
The Chern number as function of k z ; (c) The schematic plot 
of the band dispersion around the Weyl nodes in the k z =±k z 
plane, and the inset shows the chiral spin texture, (d) The 
gauge flux evaluated as Berry curvature in the (k x , k z ) plane. 



tween C — 2 and C = planes (i.e, at k z = ±fc°) are 
topologically unavoidable Weyl nodes as addressed at the 
beginning. On the other hand, however, the crossing-B 
(i.e, the closed loop in k z = plane) is just accidental 
and it is due to the presence of crystal mirror symmetry 
with respect to the k z —0 plane. The crossing-B is not as 
stable as crossing-A in the sense that it can be eliminated 
by changing the crystal symmetry. 

Using the 8 eigen states at T, we can construct an 8x8 
effective k ■ p Kane-Hamiltonian 31| . For qualitative un- 



derstanding, however, we can downfold the 8x8 Hamil- 
tonian into a simplest 2x2 model by considering the two 
basis 1 1 , | ) and \S, 



which catch the band-inversion 



nature. 



H, 



eff 



M 
Dk z k\ 



Dk z k 2 _ 
-M 



(1) 



here k± = k x ± ik y , and M = Mq — (3k is the mass term 
expanded to the second order, with parameters Mo > 
and /3 > to ensure band inversion. Since the two ba- 
sis have opposite parity, the off-diagonal element has to 
be odd in A;. In addition, the k± has to appear to con- 
serve the angular moment along z-direction. Therefore, 
to the leading order, the k z k\ is the only possible form 
for the off-diagonal element. Evaluating the eigen values 
E(k) = ±^M 2 + D 2 k 2 (k 2 + k 2 ) 2 , we get two gapless 

solutions: one is the degenerate points along the V — Z 
line with k z — ±k c z = ±^Mq//3; the other is a circle 
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FIG. 3: (Co lor online) Edge states and Fermi arcs of 
HgCr2Se4 [31]. (a) The edge states for the plane with 
fc z =0.067r. A ribbon with two edges is used, and there are 
two edge states for each edge (because C=2). The states lo- 
cated at different edges are indicated by different colors, (b) 
The calculated Fermi arcs for the (k y ,k z ) side surface. 



around the T point in the k z = plane determined from 
the equation k 2 + k 2 = Mo/ (3. They are exactly the 
band-crossings obtained from our first-principles calcula- 
tions. Due to the presence of k 2 ^ in the off-diagonal ele- 
ment [30|| , it is easy to check that Chern number C equals 
to 2 for the planes with —k z <k z <k z and k z ^ 0. The 
band dispersions near the Weyl nodes at k z = ±k z plane 
(Fig. 2(c)) are thus quadratic rather than linear, with chi- 
ral in-plane spin texture (shown in the inset of Fig. 2(c)). 
The two Weyl nodes located at ±k z have opposite chi- 
rality due to the opposite sign of mass term, and they 
form a single pair of magnetic monopoles carrying gauge 
flux in fc-space as shown in Fig. 2(d). The band-crossing 
loop in the k z — plane is not topologically unavoidable, 
however, its existence requires that all gauge flux in the 
k z —0 plane (except the loop itself) must be zero. 

This Chern semi- metal state realized in HgCr2Se4 will 
lead to novel physical consequences, which can be mea- 
sured experimentally. First, each fc z -fixed plane with 
non-zero Chern number can be regarded as a 2D Chern 
insulator, and there must be chiral edge states for such 
plane if an edge is created. The number of edge states 
is two for the case of C=2 (see Fig. 3(a)), or zero for the 
case of C=0. If the chemical potential is located within 
the gap, only the chiral edge states can contribute to the 
fermi surface, which are isolated points for each Chern 
insulating plane but nothing for the plane with C=0. 
Therefore the trajectory of such points in the (k x , k z ) sur- 
face or (k y , k z ) surface form non-closed fermi arcs, which 
can be measured by ARPES. As shown in Fig. 3(b), the 
fermi arcs end at k z =ik zl and are interrupted by the 
k z =0 plane. This is very much different from conven- 
tional metals, where the fermi surfaces must be either 
closed or interrupted by the Brillouin zone boundary. 
The possible fermi arcs has been recently discussed from 
a view point of accidental degeneracy for pyrochlore iri- 



dates [llj . Nevertheless, for the Chern semi- metal state, 



the fermi arcs should be more stable because the band- 
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FIG. 4: (Color online) Quantized Anomalous Hall effect 
in HgCr2Se4 thin film [31]. (a) The subband energy lev- 
els at r point as function of film thickness, (b) The Hall 
conductance as function of film thickness. 



crossings are topologically unavoidable. 

The QAHE, on the other hand, is an unique physical 
consequence characterizing the Chern semi-metal nature 
of HgCr2Se4, by considering its quantum- well structure. 
For 2D Chern insulators, the transverse Hall conductance 

2 

should be quantized as a xy = C— , where C is the Chern 
number. Such a quantum Hall effect without magnetic 
field has been long-pursued [12l . ll3l . ll5l ] but never achieved 
experimentally. In HgCr2Se4, considering the k z -Rxed 
planes, the Chern number C is non-zero for limited re- 
gions of k z , and this is due to the band inversion around T 
as discussed above. In the quantum well structure, how- 
ever, those low energy states around T should be further 
quantized into subbands (labeled as \H n ) and \E n ) for 
hole and electron subbands respectively), whose energy 
levels change as function of film thickness. As shown in 
Fig. 4(a), when the thickness of the film is thin enough the 
band inversion in the bulk band structure will be removed 
entirely by the finite size effect. With the increment of 
the film thickness, finite size effect is getting weaker and 
the band inversion among these subbands restores sub- 
sequently, which leads to .jumps in the Chern number or 
the Hall coefficient a xy 14]. As shown in Fig. 4(b), if the 
film is thinner than 2lA(about 2 lattice constants), the 
a xy is zero; once the film thickness is larger than the crit- 
ical thickness, we find subsequent jumps of a xy in unit of 
2e 2 /h. In fact, the strong anomalous Hall effect has been 
observed for the bulk samples of HgCr 2 Se4 [l8[ • This 
is in sharp contrast with pyrochlore iridates, where the 
anomalous Hall effect should be vanishing due to the AF 
ordering. 
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